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that are admixtures of CP-even and CP-odd configurations provide valuable in- 
formation about CKM phases and hadronic parameters. We present the general 
formalism to accomplish this task, taking also into account penguin contributions, 
and illustrate it by considering a few specific decay modes. We give particular 
emphasis to the decay B& — ► J/tpp , which can be combined with B s — ► J/ijjcj) 
to extract the B^-B^ mixing phase and - if penguin effects in the former mode 
should be sizeable - also the angle 7 of the unitarity triangle. As an interesting by- 
product, this strategy allows us to take into account also the penguin effects in the 
extraction of the B®~B® mixing phase from B s — > J ftp (j). Moreover, a discrete am- 
biguity in the extraction of the CKM angle (3 can be resolved, and valuable insights 
into S'C/(3)-breaking effects can be obtained. Other interesting applications of the 
general formalism presented in this paper, involving B& — > pp and B s c [ — * K*K* 
decays, are also briefly noted. 
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1 Introduction 



Studies of CP violation in the 5-meson system and the determination of the three angles 
a, (3 and 7 of the usual non-squashed unitarity triangle QXJ| of the Cabibbo-Kobayashi- 
Maskawa matrix (CKM matrix) are among the central targets of future B-physics 
experiments. During the recent years, several strategies were proposed to accomplish this 
task ||. In this context, also quasi-two-body modes B q — > X\X 2 of neutral 5 g -mesons 
(q G {d, s}), where both X% and X 2 carry spin and continue to decay through CP- 
conserving interactions, are of particular interest [f|, ||. In this case, the time-dependent 
angular distribution of the decay products of X\ and X 2 provides valuable information. 
For an initially, i.e. at time t — 0, present 5°-meson, it can be written as 

/(9, $, = £ 0^(t)g^(Q, $, vl/), (1) 

k 

where we have denoted the angles describing the kinematics of the decay products of 
X\ and X 2 generically by 0, $ and Note that we have to deal, in general, with an 
arbitrary number of such angles. The observables 0^ k \t) describing the time evolution of 
the angular distribution ([!]) can be expressed in terms of real or imaginary parts of certain 
bilinear combinations of decay amplitudes. In the applications discussed in this paper, 
we will focus on B q — > \X\ X 2 ]f decays, where X± and X 2 are both vector mesons, and / 
denotes a final-state configuration with CP eigenvalue rjf . It is convenient to analyse such 
modes in terms of the linear polarization amplitudes A (t), A\\(t) and A±(t) 0. Whereas 
A±(t) describes a CP-odd final-state configuration, both A (t) and A\\(t) correspond to 
CP-even final-state configurations, i.e. to the CP eigenvalues —1 and +1, respectively. 
The observables of the corresponding angular distribution are given by 

\A f (t)\ 2 with /e{0,||,±}, (2) 

as well as by the interference terms 

K{A*(t)A||(t)} and %{A}(t)A ± (t)} with / E {0, ||}. (3) 

This formalism is discussed in more detail in 0, where several explicit angular distribu- 
tions can be found and appropriate weighting functions to extract their observables in 
an efficient way from the experimental data are given. 

In the following considerations, the main role is played by neutral B q — > [XiX^j 
decays, where the "unevolved" decay amplitudes can be expressed as 



A f =Af f 



l-b f e lp fe +iuj (4) 



A f = n f Nf [l-ft/e^e - *"] , (5) 

where ui denotes a CP- violating weak phase and A/} = \Mf\e l6f . Both pj and 6f are CP- 
conserving strong phases. In this case, the observables (§) and @ allow us to probe the 
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B®-B® mixing phase cj) q and the weak phase u, as we will show in this paper. Concerning 
practical applications, u is given by one of the angles of the unitarity triangle. However, 
the observables specified in (fj) and (||) are not independent from one another and do not 
provide sufficient information to extract <fi q and uj, as well as the corresponding hadronic 
parameters, simultaneously. To this end, we have to use an additional input. 

Usually, the weak phase uj is of central interest. If we fix the mixing phase <ft q sepa- 
rately, it is possible to determine uj - and interesting hadronic quantities - as a function 
of a single hadronic parameter in a theoretically clean way. If we determine this quantity, 
for instance, by comparing B q — > X\ X2 with an S77(3)-related mode, all remaining pa- 
rameters, including uj, can be extracted. If we are willing to make more extensive use of 
flavour-symmetry arguments, it is in principle possible to determine the Bg-B® mixing 
phase 4> q as well. An example for such a strategy is given by the decay B d — > J/ifjp°, 
which can be combined with B s —>■ J/tp cf) to extract the B®-B® mixing phase (f>d = 2(3 and 
- if penguin effects in the former mode should be sizeable - also the angle 7 of the uni- 
tarity triangle. As an interesting by-product, this strategy allows us to take into account 
also the penguin effects in the extraction of the B Q S -B® mixing phase from B s — > J/ip(j), 
which is an important issue for "second-generation" 5-physics experiments at hadron 
machines. Moreover, we may resolve a discrete ambiguity in the extraction of the CKM 
angle (3, and may obtain valuable insights into S'[/(3)-breaking effects. Other interesting 
applications of the general formalism presented in this paper, involving Bd — > pp and 
B Si d — ► K*K* decays, are also briefly noted. 

As the extraction of uj with the help of these modes involves "penguin", i.e. flavour- 
changing neutral-current (FCNC) processes and relies moreover on the unitarity of the 
CKM matrix, it may well be affected by new physics. In such a case, discrepancies would 
show up with other strategies to determine this phase, for example with the theoretically 
clean extractions of 7 making use of pure "tree" decays such as B s — > DfK T . Since no 
FCNC processes contribute to the decay amplitudes of these modes, it is quite unlikely 
that they - and the extracted value of 7 - are significantly affected by new physics. 

The outline of this paper is as follows: in Section |2|, the time-dependent observables of 
the B q — > X\ X2 angular distribution are given. The strategies to extract CKM phases, as 
well as interesting hadronic parameters, with the help of these observables are discussed 
in Section ^. In Section |], we focus on the extraction of f3 and 7 from — > J /iff p° and 
B s — > J /if; cf). Further applications of the formalism developed in Sections |2| and ^ are 
discussed in Section and the conclusions are summarized in Section |6|. 

2 The Time Evolution of the Angular Distributions 

In this section, we consider the general case of a neutral quasi- two-body decay B q — ► 
[Xi Xz]f into a final-state configuration / with CP eigenvalue rjf that exhibits "unmixed" 
decay amplitudes of the same structure as those given in (§) and (H). If we use linear 
polarization states to characterize the final-state configurations as, for example, in J/J, 
we have / G {0, ||, _L}. 
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At this point a comment on the angular distribution of the CP-conjugate decay B® 
X± X 2 , which is given by 



(6) 



is in order. Since the meson content of the X\ X 2 states is the same whether these result 
from the B q or B® decays, we may use the same generic angles G, $ and \I> to describe 
the angular distribution of their decay products. Within this formalism, the effects of CP 
transformations relating B q — > [X\ X 2 ]f to B® — > [X\ X 2 ]f are taken into account by the 
CP eigenvalue r/f appearing in fl5|), and do not affect the form of g^ k \Q, <3?, Therefore 
the same functions g( k \@, are present in (Jl]) and (j^j (see also 

In view of applications to B s decays, we allow for a non-vanishing width difference 
Ar 9 = - r[ 9) between the B q mass eigenstates Bf ("heavy") and B\ ("light"). In 
contrast to the Bd case, this width difference may be sizeable in the B s system ||; it 
may allow studies of CP violation with "untagged" B s data samples, where one does not 
distinguish between initially, i.e. at time t = 0, present B® or B® mesons |J. The time 
evolution of the observables corresponding to (0) takes the following form: 



\A f (t)\ 2 

\Mt)\ 2 - 



1 
2 
1 
2 



r[ e" r L * + e" r H * + 2 e" r ^ {A f D cos(AM q t) + sin(AM 9 i)} 
R[ e- v( " h + e- r H 9) ' - 2 e~ r «* {a£ cos(AM g t) + A f M sin(AM^)} 



(7) 



where AM, = - M^ J > denotes the mass difference between the B q mass eigen- 
states, and T q = + /2. The quantities R[, R^, A^, and A^, which are not 
independent from one another and satisfy the relation 



r(«) 



(i4£) a +(Ai) 2 = i2£i^, 



(9) 



are given by 
Rl 



Wf\ 2 [(l + Vf cos0 g ) 
—2 6/ cosp/ {cosu; + 77/ cos( 
4 = |jV>i 2 [(I-77/ cos <f> q ) 



-2 6y cosp/ {cosco> — 77/ cos( 
!>/ sin pf sin a; 



+ cj)} + b) {1 + 77/ cos(0 g + 2cj)} 



+ w)} + b) {1 - 77/ cos(0 9 + 2a;)} 



2\M f \ 2 b j 



Vf l-^/l 2 s i n 0<j — 2 6/ cos pf sm((p q + to) + b 2 j sm((p q + 2cu 



Here the phase g denotes the CP-violating weak B q -B^ mixing phase: 

f 2(3 for q = d 
q I —2^7 for q — s, 



(10) 

(11) 
(12) 
(13) 

(14) 
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where 25^ ~ 0.03 is tiny in the Standard Model because of a Cabibbo suppression of 
0{\ 2 ). This phase cancels in 



S 



f 



- [rI + rQ) = (TV/I 2 (l -2b f cosp f cosu + bfj . (15) 



It is also interesting to note that there are no AM q t terms present in the "untagged" 
combination 

(16) 



whereas 



A f {t)\ 2 -\A f (t)\ 2 = 2e" r '* \A f D cos(AM q t) + A f M sin( AM q t) 



(17) 



Because of @, each of the \A f (t)\ 2 or \A f (t)\ 2 (/ G {0, terms of the B q -> X X X 2 

angular distribution provides three independent observables, which we may choose as A(^, 
A^ and Sf. 

The time evolution of the interference terms (|3|) is analogous to (0) and (|T7|) . Let us 
first give the expressions for the observables corresponding to QToT): 



R =2 



JtK(0)^,(0)} + JtK(0)A||(0)} = |JV„ 



COS A n 



|&o cos(p - A 0) ||) + fey cos(pn + A ,||)} cos^ + b b\\ cos(p - P\\ - A ,||) (18) 



where the 



%{A}(0)A ± (0)} + S{A* f (0)A ± (0)} 

b f cos(pf - A /)± ) - 6_lcos(p ± + A/,_ 



A /J = 6 f ~ S f 



sinw, 



(19) 



(20) 



denote the differences of the CP-conserving strong phases of the amplitudes Mf = e f \J\ff\ 
and Nj = e l5 f\J\fj\. On the other hand, the rate differences corresponding to ([[7]) take 
the following form: 

R{A* (t)A\\(t)} - R{jZ(tJA\\(t)} = 2e- r " t \R D cos(AM q t) + R M sm(AM q t)] (21) 



%{A}{t)A ± {t)} - %{A){t)A ± {t)} = 2e~ l " t [l f cos(AM q t) - J^sin(AM g t) 
where 



smi^j 



R D = \Af \\Af\\ \ [6osin(po-A 0l ||) + 6|| sin(pn + A ,||) 
Rm = |A/o||A/||| cos A ,|| sin0 g - {& cos(p - A ,||) + b\\ cos(p\\ + A ,||)} sm(<f) q + 



:22) 



(23) 



+ b by cos(p - P|| - A ,||) sin(0 g + 2uj) 



(24) 
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A < 

Figure 1: The amplitudes Aq, An and Aq, An in the complex plane. 



and 

// = |A//||A/j_| sin Af t ± + {bf sin(p/ — A^j.) — b± sin(pj_ + A^j.)} coscj 

- bfb±sin(p f - pj_ - A f) x) (25) 

Iyi = \Nf\\N±\ cos A/ 5 _l cos0 g - {bf cos(p/ - A/ 5 _l) + bj_ cos(p_L + A/ ; _l)} cos(4> q + u) 

+ bf b± cos(p/ — p± — Af,±) cos(0 g + 2 . (26) 

Note that / G {0, ||} in ( |T9D and fl2"2"|). The minus sign in the latter expression is due 
to the different CP eigenvalues of / G {0, ||} and / = _L. If we set "0 = || " in (|23|) and 
(p^), we get expressions taking the same form as (|T2"| ) and (|13|), which provides a nice 
cross check. The expressions given above generalize those derived in (7| in two respects: 
they take into account penguin contributions, and they allow for a sizeable value of the 
B®-B® mixing phase 4> q . In the discussion of B s — > J/xjj <f> in , it was assumed that <f) s 
is a small phase, and terms of 0(<pf) were neglected. 

Unfortunately, not all of the observables Sf, A^ and A^ are independent from those 
of the interference terms (§). This can be seen by considering two different final-state 
configurations / and /. In this case, the time-dependent angular distribution provides 
nine observables. To be definite, let us consider the case / = and / = ||. Then we 
have six observables, corresponding to Sf, A^ and A{j (/ G {0, ||}), as well as the three 
observables R, R& and Ru, which are due to the real parts in The measurement of 
Sf and Aj^ allows us to fix the magnitudes \Aq\, \A\\ \ and |Ao|, |A|||. Using in addition the 
observables R and Rd, we can determine the angle a between the unmixed amplitudes 
A and A\\, as well as the angle a between A , A\\ (see Fig. [[]). So far, the relative 
orientation of the amplitudes (Aq, An) and (Aq, Am) is not determined. However, if 
we use, in addition, the mixing-induced CP asymmetry A^, we are in a position to fix 
(fiq + ^0; where ip denotes the angle between the amplitudes A and A : 

A° M = |A o ||A o |sin(0 9 + 7/; o ). (27) 
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Since the relative orientation of the amplitudes e~ l ^ q A\\ and A\\ is also fixed this way, 
we can predict the values of the two remaining mixing-induced CP-violating observables 
A M and Ru- Consequently, only seven of the nine observables are independent from one 
another. 

It is convenient to introduce the following "normalized" observables: 



(28) 
(29) 

(30) 



R 



h = 



A f 




A f = 


A f 

V 




R 


Rd — - 


Rd 


Rm — 


Rm 








If 


ff - 


I f 


ff - 


I f 


Js f s ± 


J D — 


Js f s± 


J M — 


JSfS ± 



which have the advantage that they do not depend on the overall normalization factors 
\Aff\- The observables A^ and A M allow us to determine the hadronic parameters 6/ and 
Pf as functions of u and (p q : 



If ± Jlj - hfkf 



2b f cos p f = u f + v f b) 
2 bf sin pf — (1 — Uf cosu) + (1 — Vf cosu;) b 2 < 



'Mr 

sino; 



(31) 

(32) 
(33) 



where 



with 



and 



hf 
kf 
h 



uj + Df (1 — Uf cosu) 2 



COS LU } 



v} + D f (l-Vf 
2 — UfVf — Df (1 — Uf cosu;)(l — Vf cosu), 



u f 



v f 



rjfA^ cosu — sin(0 g + u) 

Vf A L ~ sin (0g + 2uJ ) 
rjfAljL cosu — sin(0 g + u) 



'A: 

smu 



(34) 
(35) 
(36) 

(37) 
(38) 

(39) 
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It should be emphasized that no approximations were made in order to derive these 
expressions. If we consider, in addition to A D and A^, either the observables specified 
in (29) or those given in ([30]), we obtain seven normalized observables, which depend on 



five hadronic parameters (bf, pf, bj, pj and A^ ^), as well as on the two CP-violating 
weak phases <p q and uj. However, only five of the seven observables are independent 
from one another, so that we have not sufficient observables at our disposal to extract 
these parameters simultaneously. To accomplish this task, we have to make use, for 
example, of another decay that can be related to B q — > X\ X 2 through flavour-symmetry 
arguments. On the other hand, if we use <p q and uo as an input, all hadronic parameters 
describing the decay B q — > X\ X 2 can be extracted without any additional assumption, 
thereby providing valuable insights into hadronic physics and a very fertile testing ground 
for model calculations of B q — > X\ X 2 . The measurement of the angular distributions 
discussed in this paper requires high statistics and can probably only be performed at 
"second-generation" S-physics experiments at hadron machines, such as LHCb or BTeV, 
where also decays of 5 s -mesons can be studied. Since several promising strategies to 
extract the weak phases (fi q and uj at such experiments were already proposed (see, for 
example, |§), it may indeed be an interesting alternative to use measurements of angular 
distributions not to extract CKM phases, but to explore hadronic physics. 

In practical applications, the parameters bf typically measure the ratio of "penguin" 



to "tree" contributions. Applying the Bander-Silverman-Soni mechanism [10], and fol 



lowing the formalism developed in |TT|, 12], which makes use - among other things - of 



the "factorization hypothesis" , we obtain for various classes of B decays 

b f = b, p f = p V/e{0,||,_L}. (40) 

The main reason for these relations is that the form factors, which depend on the final- 
state configuration /, cancel in the ratios bf of "penguin" to "tree" contributions. Al- 
though non-factorizable contributions are expected to play an important role, thereby 
affecting (fH]), it is interesting to investigate the implications of these relations on the 
observables of the angular distributions in more detail. If we introduce 

i _ 2 6sinpsina; 

D 1 — 2 b cosp cos a; + b 2 

- sin (f) q - 2 b cos p sm((j) q + u) + b 2 sin(0 g + 2 uj) 

Am ~ 1 ol m ' 

1 — lb cos p cos uj + b z 



we obtain 



A f B = A D , A f M = r) f A M , (43) 
-R = cosA 0i ||, i?D = Ad cos Aoj, -Rm = Am cos Aoj, (44) 
J / = sinA /i± fl = A D sin A ft± , (45) 



ff 



cos Ajj_ 



cos 4> q — 2b cos p cos(0 g + uj) + b 2 cos(0q + 2 uj) 



1 — 2 b cos p cos uj + b 2 



(46) 
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where 

, 2 



(A D ) +(A m ) +I 2 = 1. (47) 



These relations provide an interesting test of whether (fiOf) is realized in the decay B q 
Xi X 2 . Note that does not - in contrast to (pE5|) - vanish for trivial values of A^j_. 



3 Extracting CKM Phases and Hadronic Parameters 

Let us now focus on the extraction of CKM phases from the observables of the B q — > 
Xi X2 angular distribution. As we have already noted, to this end, we have to employ an 
additional input, since we have only five independent normalized observables at our dis- 
posal, which depend on seven "unknowns" . Although it would be desirable to determine 
(f) q and u> simultaneously, usually only the CKM phase us is of central interest. 

The B®-B® mixing phase </> s = —25^ = 2 axg(V t * V t b) is negligibly small in the Standard 
Model. It can be probed - and in principle even determined - with the help of the decay 
B s — > J/tj) (f) (see, for example, 0]). Large CP- violating effects in this decay would signal 
that 2^7 is not tiny, and would be a strong indication for new-physics contributions to 
B° s -B® mixing. On the other hand, the B^-B® mixing phase <pd = 2/5 can be fixed 
in a reliable way through the "gold-plated" mode Bd — > J/tpKs Hl3| . Strictly speaking, 



mixing- induced CP violation in Bd — > J /ijj K§ probes sin(2 + where 4>k is related to 
the weak K°-K° mixing phase and is negligibly small in the Standard Model. Because 
of the small value of the CP-violating parameter ek of the neutral kaon system, <pK 
can only be affected by very contrived models of new physics [Q. A measurement of 
mixing- induced CP violation in Bd — > J/ipK^ allows us to fix <pd = 2/5 only up to a 
twofold ambiguity. Several strategies to resolve this ambiguity were proposed in the 
literature flI5| , which should be feasible for "second-generation" 5-physics experiments. 
As we will see in the following section, also the decay Bd — > J/ip P°, i n combination with 
B s — > J/ip </>, allows us to accomplish this task. 

If we use <j) q thus determined as an input and consider, in addition to A D and A{j, 
either the observables specified in fl29|) or those given in (|3~0D, we can determine u as 
a function of a single hadronic parameter. Let us, for the moment, focus on the latter 
case, i.e. on the observables A D , A^, A^, A^ and //, J D , for a given final-state 
configuration / £ {0, ||}. Since |JV/| and \N±_\ cancel in these quantities, they depend 
only on the hadronic parameters bj, pf, p±, Af t ±, as well as on the weak phases 
u and (f) q . Consequently, we have seven observables at our disposal, which depend on 
seven "unknowns" . However, only five of the seven observables are independent from one 
another, as we have discussed in the previous section. If we use <p q as an input, we can, 
for instance, obtain uj and bf, pf, b±, p± as functions of the strong phase difference A^j, 
in a theoretically clean way. Although the following discussion deals with Ajj_, we can 
also replace this quantity by another hadronic parameter of our choice. If we fix A^j_, for 
example, by comparing B q — > X\ X2 with an 5 , t/(3)-related mode, all parameters can be 
extracted. Using in addition the observables Sf, we can also determine the normalization 
factors \Nf\. Comparing them with those of the S'L r (3)-related mode used to fix Ayy., 
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we can obtain valuable insights into S77(3)-breaking corrections. The observables that 
we have not used so far can be used to resolve discrete ambiguities, arising typically in 
the extraction of these parameters. 

Let us now give the formulae to implement this approach in a mathematical way. The 
general expression for the observable If (see fl25D and (|3~0"D) leads to the equation 



Afsm Af.i + Bf cos A 



c 



(46 



where 
Af 

Bf 

with 



N, 
1 

w f 



1 — (bf cos pf + b± cos p±) cosu + bfb± (cos pf cos p± + sinpj sin p_|_) (49) 



(bf sin pf — b± sin p±) cos to — 6/6^ (sin pf cos p_i_ — cos pf sin p_|_) 



The solution of 



1 — 26/ cos pf cos a; + 6jJ (1 — 2 6j_ cos p± costo + 6j 
is straightforward, and is given as follows: 



(50) 
(51) 

(52) 



sin A 



A f Cf± 



/,-L 



A 2 



R2 
^7 



ci Kb 



Aj + Bj 



cos A 



C/ - A f sin A 



/,-L 
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/ 



(53) 



If we insert 6/ and p/, determined as functions of and (fi q with the help of ( pT]) -([53]) , 
into the expressions given above, we can - for a given value of cj) q - determine Aj j_ as a 
function of to. It should be emphasized that the relation between A^j_, to and <p q obtained 
this way is valid exactly. Using or instead of If would lead to the same relation, 
since these observables are not independent from If. 

Alternatively, we may use the observables (|29]) instead of (|3~0"D . The general expression 
for R (see (|8|) and ( p9|) ) implies an equation similar to p8|) , where Af,Bf and C/ have 
to be replaced through 



.4 
C 



iV 

1 

iV 



sin p|| — 6 sin p ) cos a; + 6 ^|| (sin Po cos P|| — cos Po sm P| 



(&o cos Po + &|| cos P|| ) cos w + 6q^|| (cos po cos p\\ + sin p sin p| 



(54) 

(55) 
(56) 



where 



N = y (l — 2 6 cos po cos c<j + &§) (l — 2 6|| cos py cos a; + 6|) . 



(57) 



Obviously, the most efficient strategy of combining the observables provided by the B q 
Xi X 2 angular distribution depends on their actually measured values. 
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If we are willing to make more extensive use of flavour-symmetry arguments than 
just to fix the strong phase difference A?y, it is in principle possible to determine also 
the Bq-B® mixing phase <f> q . In the following section, we will have a closer look at the 
decay Bd — > J/ipp°, which can be related to B s — > J/ip(f> through SU(3) arguments 
and a certain dynamical assumption concerning "exchange" and "penguin annihilation" 
topologies. However, before we turn to these modes, which allow the simultaneous ex- 
traction of (ftd = 2/9 and 7, let us first give two useful expressions for the observables R 
and If. Since the parameters bf measure typically the importance of "penguin" topolo- 
gies in comparison with current-current contributions, they may not be too large. If we 
eliminate the hadronic parameters bf and pf in R and If with the help of the observables 
Ajy and and keep only the leading-order terms in bf, we obtain 

« K cosA ,„-I(i»-4)^ (58) 

^sinA^ + i^-ii)^, (59) 

allowing us to determine u if the strong phase differences A 0j || or A^j_ are known. Inter- 
estingly, the leading-order expressions (|58|) and ( |59D do not depend on the B®-B® mixing 

phase (f) q . A possible disadvantage of R is that u enters in combination with sinAo,||. 
Since A 0) || is a difference of CP-conserving strong phases (see (|20|) ), it may be small, 
thereby weakening the sensitivity of these observables on u. The situation concerning 
this point is very different in the case of the observables If, which allow us to determine 
u even in the case of A ft± G {0°, 180°}. 



4 Extracting (3 and 7 from Bd — > J/iJj p° and 

If we combine the observables describing the time-dependent angular distribution of the 
decay B d -> J/V>[-> tt+tt-] with those of B s -> J/^[-> l + l~](f>[-^ K+K'}, 

we may extract the B^-B® mixing phase <pd = 2(3 and the angle 7 of the unitarity 
triangle. The B^, — > J /if} p° angular distribution can be obtained straightforwardly from 
the B s — > J /if) (p case, which has been discussed in detail in [0] , by performing appropriate 
replacements of kinematical variables. 

The decay B® —>■ J /if) p° originates from b — > ccd quark-level transitions; the structure 
of its decay amplitude is completely analogous to the one of B® — > J/ip K$ (see [0). For 
a given final-state configuration / with CP eigenvalue rjf, we have 

A(B° d - [J ft p%) = A» [Aif + AW] + \WA$J + X^Afl , (60) 

where A^f is due to current-current contributions, and the amplitudes A^ describe 
penguin topologies with internal q quarks (q G {u,c,t}). These penguin amplitudes take 
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into account both QCD and electroweak penguin contributions, and 



= V qd V; b 



(61) 



are the usual CKM factors. Employing the unitarity of the CKM matrix and the Wolfen- 

tain 

(62) 

Aj ^XrA .4;.;:" + .4;;;:,'/ . (63) 



stein parametrization [0, generalized to include non- leading terms in A [18], we obtain 

A(B° d - [J/^p ]/) = -XA f [l - afjOf^ 

where 



with 4 c 5^41/-41/,and 



a f e»f = R* (l-y 



A( ut )f 

pen 



^lcc "T ^Ipen 



(64) 



The quantity Al^f is defined in analogy to A^f , and the relevant CKM factors are 
given as follows: 



A= \V US \ = 0.22, A= ^\V cb \ = 0.81 ±0.06, R b = j 



ub 



V, 



rb 



0.41 ±0.07. (65) 



It should be emphasized that (|62f) is a completely general parametrization of the B® — > 
J/V p° decay amplitude within the Standard Model, relying only on the unitarity of the 
CKM matrix. In particular, this expression takes also into account final-state-interaction 
effects, which can be considered as long-distance penguin topologies with internal up- 
and charm-quark exchanges | i~9fl . Comparing ( |6^) with @, we observe that 



J\ff = —\Af, bf = af, pj = 6j, to = 7. (66) 
Let us now turn to B Q S — > J/ip<f). Using the same notation as in (j62|), we have 



A(5 S O ^[J/V>0]/) 



1 i ' iO't ii 

1 + eaj-e fe 7 



(67) 



where Alt and a'^e f take the same form as fl63|) and (@), respectively, and 



e = 



1 - A 2 



(6c 



The primes remind us that we are dealing with a6->s transition. Consequently, if we 
compare (|67D with (§]), we obtain 



A^ 



l — I b f = ea' f , P f = 0' f + 180°, u = 7. 



(69) 
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The B a — > J /ip (j) and Bd ^ J/ip p° observables can be related to each other through 



\A'\ = V2\A 



f\ 



(70) 



7,/ " Hf ( 71 ) 

a' f = a f , B) = B h (72) 

where the factor of \/2 is due to the p° wave function. These relations rely both on the 
SU(3) flavour symmetry of strong interactions and on the neglect of certain "exchange" 
and "penguin annihilation" topologies. Although such topologies, which can be probed, 
for example, through B s — > p + p~, D* + D*~ decays, are usually expected to play a very 
minor role, they may in principle be enhanced through final-state-interaction effects 
For the following considerations, it is useful to introduce the quantities 



1 Sf 1 — 2 cif cos8f cos7 + a 2 

1 \\Af\ ) 5} ~ 1 + 2ea' f cos6' f cos7 + e 2 a 



12 ' 
/ 



(73) 



which can be fixed through the "untagged" Bd — > J/tpp° and B s — > J/i/j(f> observables 
with the help of ([7D]). Consequently, each of the linear polarization states / G {0, ||, J_} 
provides the following three observables: 



Hfi ^E)' ^m- 



(74) 



Applying ( [72] ) to flTB|), these observables depend only on the hadronic parameters aj and 
9f, as well as on the B^-B® mixing phase 0^ = 2(3 and the angle 7 of the unitarity 
triangle. If we choose two different linear polarization states, the observables ([74]) allow 
us to determine the corresponding hadronic parameters and /3 and 7 simultaneously. 

This approach can be implemented in a mathematical way as follows: if we consider 
a given final-state configuration / and combine the observables Hf and A^, which do not 
depend on 0^, with each other, we can determine a,f and 9f as functions of 7: 



a f 

2a/ cos9f 
2 a/ sin 9f - 



Pf ± yPf - if 



1-H f + (1- e 2 H f ) a) 



(1 + tHf) cos 7 
[l + e)(l + eaj)H/ 



[1 + eH 



f) 



sin 7 



where 



Pf 



2 (1 + eHff cos 2 7 - (1 - H f ) (1 - e 2 H f ) sin 2 7 - e f E f 



(If 



{1 -e 2 H f Y sin 2 ~f + e 2 E f 

(1 -Hf) 2 sin 2 7 + E f 
'l-e 2 H f ) 2 sm 2 j + e 2 E/' 



(75) 
(76) 
(77) 

(78) 
(79) 
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with 

cos 7 . (80) 

These expressions allow us to eliminate the hadronic parameters af and Of in the mixing- 
induced CP asymmetry A^, thereby fixing a contour in the "/-(fid plane, which is related to 

Afsm<f) d + B f cos <p d = C f , (81) 

with 

Af = 1 — 2 df cos Of COS7 + aj cos 27 (82) 

5/ = — 2 a/ cos Of sin 7 + sin 27 (83) 

C/ = (l - 2a/ cos6»/cos7 + af) (7]fA f M ). (84) 



The solution of ([H]) has already been given in (53). If we consider two different final-state 
configuratons / and /, we obtain two different contours in the "/-(fid plane; their intersec- 
tion allows us to determine both 7 and (fid = 2/?. Using, in addition, the observables 029D 
or ( |30|) - depending on which final-state configurations / and / we consider - we may 
resolve discrete ambiguities, arising typically in the extraction of (fid and 7. 

Because of the strong suppression of a'f through e = 0.05 in ([73]), this approach is 
essentially unaffected by possible corrections to ([72]), and relies predominantly on the 
relation (|70|). If we insert the values of (fi d and 7 thus determined into the expressions 
for the observables of the third linear polarization state /', which has not been used so 
far, its hadronic parameters \Af\, af and Of can also be determined. Comparing \Af\ 
with the B s — > J/ifi(fi parameter \A'f\, we can obtain valuable insights into the validity 
of (f70f ). Moreover, several other interesting cross checks can be performed with the 
many observables of the angular distributions. Because of our poor understanding of the 
hadronization dynamics of non-leptonic B decays, only the "factorization" approximation 
can be used for the time being to estimate factorizable 5*?7(3)-breaking corrections to 
( ffOj) . Explicit expressions for the B s — > J /if) (fi observables can be found in [^] , and 
SU (3)-breaking effects in the corresponding form factors were studied in |2T[ . However, 



also non-factorizable effects are expected to play an important role, and experimental 
insights into these issues would be very helpful to find a better theoretical description. 

The simultaneous extraction of (fid and 7 discussed above works only if the hadronic 
parameters a/ and Of are sufficiently different from each other for two different final- 
state configurations /. If, for example, (|40| ) should apply to Bd — > J /if) p° - which seems 
to be quite unlikely - the B^-B® mixing phase has to be fixed separately in order to 
determine 7. In this case, each linear polarization state / e {0, ||, _L} provides a strategy 
to extract 7 that is completely analogous to the one proposed in which makes use 
of B s (d) — > J I if) K$ decays. If we combine Hf with A^, we obtain 



a f 



Hf-l+u f (l + eH f )cosj ^ 



\ l-v f (l + eH f ) cos 7 - e 2 H f 
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The intersection of the contours in the 7— a/ plane described by this expression with those 
related to (|5ID allows us to determine 7 and af. 

If we use <pd as an input in order to extract 7 from Bd — > J/ifip , it is, however, 
more favourable to follow the approach discussed in the previous section, i.e. to use (|5TS|), 
and to fix Aj j_ (or A 0) |i) through the the B s — > J/ip(fi observables with the help of (|7l| ) 
and (|T2]). Using in addition the observables involving the third linear polarization state 
/' that we have not employed so far, we can also fix its hadronic parameters ay and 
8f/, as well as the strong phase difference Afij. Comparing Af/j with its B s — > J/if)(f> 
counterpart A'j, p we may obtain valuable insights into possible corrections to (|7l|). 

As an interesting by-product, this strategy allows us to take into account also the pen- 
guin effects in the extraction of the B®-B® mixing phase (fi s from B s — > J/ if) (fi. Although 
the penguin contributions are strongly suppressed in this mode because of the tiny pa- 
rameter e = 0.05 (see QB7p), they may well lead to uncertainties of the extracted value of 
(fi s at the level of 10%, since (fi s = (9(0.03) within the Standard Model. A measurement of 
(fi s = —2X 2 f) would allow us to determine the Wolfenstein parameter 77 |T7J, thereby fixing 
the height of the unitarity triangle. Since the decay B s — ► J/if)(fi is very accessible at 
"second-generation" 5-physics experiments performed at hadron machines, for instance 
at LHCb, it is an important issue to think about the hadronic uncertainties affecting the 
determination of (fi s from the corresponding angular distribution. The approach discussed 
above allows us to control these uncertainties with the help of Bd — > J/if> p°. 

The experimental feasibility of the determination of 7 from the Bd — > J /if) p° angular 
distribution depends strongly on the "penguin parameters" aj. It is very difficult to esti- 
mate these quantities theoretically. In contrast to the "usual" QCD penguin topologies, 
the QCD penguins contributing to B d — > J /if) p° require a colour-singlet exchange, i.e. are 
"Zweig-suppressed" . Such a comment does not apply to the electroweak penguins, which 
contribute in "colour-allowed" form. The current-current amplitude A)£)f originates from 
"colour-suppressed" topologies, and the ratio A^f / A^}f + A^/ , which governs a/, 
may be sizeable. It would be very important to have a better theoretical understanding 
of the quantities a / e l9f . However, such analyses are far beyond the scope of this paper, 
and are left for further studies. 

If the parameters aj should all be very small, which would be indicated by A{y = 
Rd — It> — 0) we could still determine the B^-B® mixing phase from the observables 
Ay[ = rjf sin^rf. If we use, in addition, 1^ = cos A^j_ cos (fid and fix cos Af t ± through the 
corresponding B s — > J/ip (fi observable, cos (fi d can be determined as well. Consequently, 
the Bj-B® mixing phase (fid can be fixed unambiguously this way, thereby resolving a 
twofold ambiguity, which arises in the extraction of (fid from Bd — > J/ip K$. This mode 
probes only sin (fi d . Since (fid = 2/3, we are left with a twofold ambiguity for j3 G [0°, 360°]. 
If we assume that (3 G [0°, 180°], as implied by the measured value of Sk, we can fix j3 
unambiguously. For alternative methods to deal with ambiguities of this kind, see [|1^ . 

Before we turn to Bd — > pp and B s — > K* K* decays, let us note that the approach 
presented in this section can also be applied to the angular distributions of the decay 
products of B s{d) -> J/VH l + n K*[-+ ir°K s ] and B d{s) -> D*+ s) D*^. For the B s{d) -> 
J /if) Ks and B d ( s ) — > ^d(s)^d(s) variants of these strategies, see [16|. 
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5 Further Applications 



In this section, we discuss further applications of the general strategies presented in Sec- 
tion |3|. All of the methods discussed below have counterparts using Bd,s decays into two 
pseudoscalar mesons. If we replace the pseudoscalars by higher resonances, for example, 
by vector mesons, as in the following discussion, the angular distributions of their de- 
cay products provide interesting alternative ways to extract CKM phases and hadronic 
parameters, going beyond the B d s — > PP strategies. Because of the many observables 
provided by the angular distributions, we can, moreover, perform many interesting cross 
checks, for example, of certain flavour-symmetry relations. 



5.1 The Decays B d 



p+p- 



and B< 



K*+K* 



The decay B® — > p + p originates from b — > uud quark-level processes. Using the same 
notation as in (1621), we have 



A(B° d - [p + p-] f ) 



— C f e* 7 



d f e i@ fe- h 



(86) 



where 



and 



C f = \ 3 AR b 



d 



f 



A( u )f _|_ A( ut )f 
cc ' pen 



pen 



J7) 



\ 2 /2)R b 



^lec 



I I(«t)/ 

~t~ ^Ipen 



In order to distinguish the B® — > p + p~ amplitudes from the B d —>■ J/ip p° case dis- 
cussed in the previous section, we have introduced the tildes. The phase structure of the 
B® — > p + p~ decay amplitude given in (86), which is an exact parametrization within the 
Standard Model, is completely analogous to the one for the B® — > 7r + 7r~ amplitude given 



in [221, where a more detailed discussion can be found. 



The expressions for the observables describing the time evolution of the angular dis- 
tribution of the decay products of B d — > p + [-^ ty + ty ] p~[-^ 7r~7r°] can be obtained 
straightforwardly from the formulae given in Section |^, by performing the following sub- 
stitutions: 

■A/^fi-^W bf = df, Pf = Qf, w=-7- (89) 



Because of the factor of e n in front of the square brackets on the right-hand side of 
we have to deal with a small complication. Since the observables are governed by 



1 - d f e ie fe +i ^ 



d f e^fe'^ 



we have to do the following replacement, in addition: 

<\>d -> <Pd + 27 . 



(90) 



(91) 
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Let us now turn to the decay B° s — > K* + K*~, which is due to b — > uus quark- 
level transitions. For a given final-state configuration / of the K* + K*~~ pair, its decay 
amplitude can be parametrized as follows: 



A{B° S -> [K* + K* 



XC' f e 



'7 



1 + 



1-A 2 
A 2 



7/ ie' f - 

d f e fe 
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(92) 



where C'j and d'^e 10 ? take the same form as flB7|) and (0), respectively, and the primes 
have been introduced to remind us that we are dealing with afc^s mode. The phase 
structure of ( |9"2] ) is completely analogous to the B® — > K + K~ decay amplitude |2"2]| . The 
observables of the time-dependent B s — > K* + [— > irK] K*~ [— > WK) angular distribution 
can be obtained straightforwardly from the formulae given in Section |2] by simply using 
the replacements: 



Afj = AC' 



1-A 2 



/• 



d' f , Pf = & f + 180°, 



UJ 



-j. 



(93) 



Moreover, we have to perform the substitution 



<P S - 4> s + 2 7 (94) 

because of the factor of e n in front of the square brackets in ( p2|) . 

Explicit expressions for the Bd — > p + p~ and 5 S — > K* + K*~ angular distributions in 
terms of helicity amplitudes can be found in [ I2"| . Since B^ —>■ p + p~ and — > K* + K*~ 
are related to each other by interchanging all down and strange quarks, the £/-spin flavour 
symmetry of strong interactions implies 



\C>, 



/i 



\c 



/I' 



d' f = df, & f = 



.f 



as well as 



A': 



A, 



(95) 



(96) 



7,/ ~ 

In contrast to (|70|)-(|T2|), these relations do not rely on any dynamical assumption - just 
on the [/-spin flavour symmetry. They can be used to combine the B^ — > p + p~ and 
B s — > K* + K*~ observables with each other, thereby allowing the extraction of the CKM 
angle 7 and of the B^ s -B^ s mixing phases (pd = 2/3 and <p s = —28^. In contrast to 
the Bd — > 7r + 7r~, B s — > K + K~ variant of this approach proposed in ||22|| , both mixing 



phases and the CKM angle 7 can in principle be determined simultaneously. However, 
for the extraction of 7, it is more favourable to fix (pd and S separately. Then we are in 
a position to determine two contours in the 7-A j j and 7~A'^ planes in a theoretically 
clean way with the help of (|53|) . Using now the [/-spin relation (p6|) , 7 and all hadronic 
parameters describing the decays Bd — > p + p~ and B s K* + K*~ can be determined. 
As we have already noted, the hadronic parameters provide a very fertile testing ground 
for model calculations of the decays B d — > p + p~ and B s — > K* + K*~ . In particular, 
the penguin parameters dfe ie f and d' f e t0 f would be very interesting; comparing their 
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values with each other, we could obtain valuable insights into [/-spin-breaking corrections. 
Moreover, there is one strong phase difference A//j left, which can be compared with its 
[/-spin counterpart A'f, j. If we should find a small difference between these phases, it 
would be quite convincing to assume that our [/-spin input (^) is also not affected by 
large corrections. 

Let us finally note that there is another interesting way to parametrize the — > p + p~ 
decay amplitudes (see also p3|). If we eliminate through the unitarity of the CKM 

Ad) 



matrix - instead of \ t , as done in (|8q) - we obtain 



A(B° d - [p+p- 




^cc ~ ^pen 



I + 77 e la f'e 



(97) 



where 



with 



r f e 



1 - X 2 /2)R b 



pen 



!(«)/ 

^lcc 



T ^Ipen 



Rt, = T- 



V* 



td 



0(1). 



2/3 and /3 + 7 



Taking into account that we have cj)d 
Model, we arrive at 

b f = r /> Pi = ® i 1 



(98) 



(99) 



180° — a within the Standard 



to 



a, 



(100) 



and at the "effective" mixing phase = ^+27 = —2a. Consequently, using the strategy 
presented in Section |3], the 5^ — ► p + p~ angular distribution allows us to probe also the 
combination a = 180°— (3— 7 directly, i.e. to determine a as a function of a CP-conserving 
strong phase difference Aj^. (see also ([58]) and Q59|)). Needless to note that the decay 
5^ — ► p°p° may also be interesting in this respect. Since the normalization factors Mf of 
the parametrization (|9~7"D are proportional to 



A( u )S _|_ A( uc )f 
cc ~ pen ) 



(101) 



which is governed by "colour-allowed tree-diagram-like" topologies, it may well be that 
A 



0. This relation would allow us to extract a, as well as the — > p + p~ hadronic 
parameters, which include also another strong phase difference Apj, providing an im- 
portant cross check. 



5.2 The Decays B d -> K*° K*° and B s K* Q K*° 

The decays — > and S° — > are pure "penguin" modes, originating 

from b — > dss and 6 — > sdd quark-level transitions, respectively. They do not receive 
contributions from current-current operators at the "tree" level, and can be parametrized 
within the Standard Model in complete analogy to ( |62]) and fl67|). We have just to set the 
current-current amplitudes equal to zero in these expressions. The decays B® — > K*° K*° 
and .8° — > K*° K*° are related to each other by interchanging all down and strange 
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quarks, i.e. through the [/-spin flavour symmetry of strong interactions, and the strategies 
to probe 7 and the B% s -B® s mixing phases are analogous to those discussed in Section |. 
Since the Bd tS — * K*° K*° decays are pure "penguin" modes, they represent a particularly 
sensitive probe for new physics. 

An interesting alternative to parametrize the B% 
within the Standard Model is given as follows: 



[K*° K*°]f decay amplitudes 



A(B [ > -> [K*° K*°]f 



9/e 



i<Pf p if3 



where 



9je 



1 A( cu )f 



(102) 



(103) 



may well be sizeable due to the presence of final-state-interaction effects EM. Conse- 
quently, we have 

b f = 9f, Pf = <Pf, w = (3. (104) 

Because of the factor of e~ l @ in front of the square brackets on the right-hand side of 
( |102| ), the "effective" mixing phase is given by = <pd — 2(3. Consequently, the strategy 
presented in Section [| allows us to probe the CP-violating weak phase (3 of the CKM 
element V t d = \V t d\e~ 1 ^ ■ Within the Standard Model, we have <fi = <pd — 2(3 = 0. However, 
this relation may well be affected by new physics, and represents a powerful test of the 
Standard-Model description of CP violation (for a recent discussion, see PJ). Therefore it 
would be very important to determine this combination of CKM phases experimentally. 
The observables of the B d — > K*°[-^ n~K + ] K*°{-^ n + K~] angular distribution may 
provide an important step towards this goal. 



6 Conclusions 

The angular distributions of certain quasi-two-body modes B^^ XiX 2 , where both 
X\ and X 2 carry spin and continue to decay through CP-conserving interactions, provide 
valuable information about CKM phases and hadronic parameters. We have presented 
the general formalism to accomplish this task, taking into account also penguin contri- 
butions, and have illustrated it by having a closer look at a few specific decay modes. In 
comparison with strategies using non-leptonic B d s decays into two pseudoscalar mesons, 
an important advantage of the angular distributions is that they provide much more 
information, thereby allowing various interesting cross checks, for instance, of certain 
flavour-symmetry relations. Moreover, they provide a very fertile testing ground for 
model calculations of the B^ s X\ X 2 modes. 

We have pointed out that the decay B d — > J ftp p° can be combined with B s — ► J/ip <f> 
to extract the B^-B® mixing phase <pd = 2(3 and - if penguin effects in the former 
mode should be sizeable - also the angle 7 of the unitarity triangle. As an interesting 
by-product, this strategy allows us to take into account also the penguin effects in the 
extraction of the B° s -B® mixing phase from B s — ► J /ip (f). If penguin effects should be very 
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small in Bd — > J/ip p°, (fid could still be determined and it would even be possible to resolve 
a twofold ambiguity, arising in the extraction of this CKM phase from Bd — > J/ifiKg. 
Other interesting applications, involving Bd — > pp and B s d — > K*K* decays, were also 
noted. Within the Standard Model, these modes are expected to exhibit branching ratios 
at the 10~ 5 level; also the one for B d — > K* ^* may well be enhanced, from its "short- 
distance" expectation of 0(1O~ 6 ) to this level, by final-state-interaction effects. 

Since the formalism presented in this paper is very general, it can of course be ap- 
plied to many other decays. Detailed studies are required to explore which channels are 
most promising from an experimental point of view. Although the Bd modes listed above 
may already be accessible at the asymmetric e + -e~ B-factories operating at the T(4S) 
resonance, which will start taking data very soon, the strategies presented in this paper 
appear to be particularly interesting for "second-generation" experiments at hadron ma- 
chines, such as LHCb or BTeV, where also the very powerful physics potential of the B s 
system can be exploited. 
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